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Xinyue Diao,':* Fenhai Guan,! Yijie Wang,! Yuhao Qin,! Zhi Qin,! Dong Guo,! Qianghua Wu,! Dawei 
Si,! Xuan Zhao,! Sheng Xiao,! Yaopeng Zhang,! Xianglun Wei,” Haichuan Zou,* Herun Yang,” Peng 
Ma,” Rongjiang Hu,” Limin Duan,” Artur Dobrowolski,* Krzysztof Pomorski,* and Zhigang Xiao! 


! Department of Physics, Tsinghua University, Beijing 100084, China 
? Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, China 
School of Nuclear Science and Technology, Lanzhou University, Lanzhou 730000, China 
“Uniwersytet Marii Curie Sktodowskiej, Katedra Fizyki Teoretycznej, Lublin 20031, Poland 


We report the reconstruction method of the fast fission events in 25 MeV/u ®°Kr +?°8Pb reactions at the 
Compact Spectrometer for Heavy IoN Experiment (CSHINE). The fission fragments are measured by three 
large-area parallel plate avalanche counters, which can deliver the position and the arrival timing information 
of the fragments. The start timing information is given by the radio frequency of the cyclotron. Using the 
velocities of the two fission fragments, the fission events are reconstructed. The broadening of both the velocity 
distribution and the azimuthal difference of the fission fragments decrease with the folding angle, in accordance 
with the picture that fast fission occurs. The anisotropic angular distribution of the fission axis also reveals 


consistently the dynamic feature the fission events. 
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I. INTRODUCTION 


One of the purposes of studying the heavy ion reactions 
(HIR) is to infer the properties of nuclear equation of state 
(EOS), which is an essential input in modeling and computing 
the evolution and properties of neutron stars and their merging 
[1-3]. The isovector sector of the nuclear EOS, namely the 
density behavior of the symmetry energy Esym(p), has been 
a long-standing open question in nuclear physics. At Fermi 
energies, some observables have been identified to constrain 
Esym(p) near normal density, including isospin diffusion, 
dipole polarizability and particle emissions etc [4-10]. Very 
recently, the PREX II experiment has reported the result of the 
neutron skin thickness, which yields a stiff E;,,,(p) in ten- 
tion with previously existing constraints [11, 12]. At supra- 
saturation densities, many experiments have been progressed 
greatly to measure the charged pion ratios or collective flow in 
heavy ion collisions to probe E;,,,(p) [13, 14]. An External- 
target Experiment on HIRFL-CSR (CEE) is under construc- 
tion for the studies on this direction [15, 16]. 

Nuclear fission is a large-amplitude collective motion 
mode involving up to hundreds of nucleons. Recently the 
studies on nuclear fission have been revived for its signifi- 
cance in both nuclear physics and astrophysics. In the stel- 
lar environment, the abundance of the nuclide in A œ~ 160 
region is significantly influenced by the recycling of the fis- 
sion products [17—19]. Theoretically, statistical fission has 
been described well by microscopical theories and various 
phenomenological approaches [20-32]. When the excita- 
tion energy or the angular momentum becomes much high, 
as achieved in HIRs well above the Coulomb barrier, the 
fission barrier tends to vanish. As a consequence, the fis- 
sion time scale becomes shorter by a factor of 10 to 100 
and the variance of the mass asymmetry increases signifi- 
cantly (the mass asymmetry 7 of the two fragments can be 
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larger than 0.6), compared to the statistical fission [33-37]. 
In this case, the dynamic feature of the fission process is of 
significance and the transport models have been successfully 
applied to describe the fast fission process [38—46].Time- 
dependent Hartree-Fock theory can also give an excellent de- 
scription of fast fission starting with large deformation [47]. 
Experimentally, fast fission, usually termed as dynamic fis- 
sion because of its vanishing fission barrier and short time 
scale, has been investigated in various systems in the last three 
decades [48—54]. 


The topic of fast fission with the simultaneous emission of 
the particles deserves further investigations because the fis- 
sioning system provides an appropriate laboratory to probe 
Esym(P). The connection between the studies of fast fission 
and Esym(p) has been recently established by the simula- 
tions with improved quantum molecular dynamics (ImQMD) 
[55, 56]. It is suggested that fast fission process follow- 
ing HIR carries the effect of Esym(e) and provides sensi- 
tive probes, because of the formation of the low-density and 
neutron-rich neck and of larger surface of two fragments com- 
pared to the non-fission process [56, 57]. The isospin content 
of the light particles emitted from the fast fission events has 
been used to probe Esym(p) experimentally [10]. 


In order to conduct the experimental studies of the fast fis- 
sion and the coincident emission of light charged particles 
(LCP) and intermediate mass fragments (IMFs), a compact 
spectrometer for heavy ion experiment (CSHINE) has been 
built [58, 59]. While the LCPs and IMFs are measured using 
the silicon strip detector telescopes (SSDT) [60, 61], the fis- 
sion fragments are measured by the parallel plate avalanche 
counters (PPACs) [62]. 


In this paper, we present the measurement of fast fission 
in 8°Kr +7°8Pb reactions with CSHINE in its second phase. 
After a brief introduction to the phase-II setup of CSHINE in 
section II, the reconstruction of velocity of the FF will be in- 
troduced in Section III and the dynamic feature of fast fission 
are presented in Section IV. Section V is the summary. 


Il. CSHINE DETECTOR SYSTEM AND THE 


EXPERIMENTAL SETUP 


72 
73 


7 The beam experiment was performed at the Radioactive 
7s Ion Beam Line I (RIBLL1) on the Heavy Ion Research Fa- 
z cility in Lanzhou (HIRFL), China. The ?95Pb target of 1 
7 mg/cm? areal density was bombarded with 25 MeV/u °°Kr 
7s beam. The charged reaction products are measured using 
77 CSHINE, which are installed in a large scattering cham- 
so ber located at the final focal plane of RIBLL1. In the cur- 
rent experiment, three PPACs are installed for fission frag- 
ment measurements to reconstruct the reaction geometry. Be- 
sides, four SSDTs are installed covering the polar angle range 
10° < flab < 60°. The SSDT is a three-layer detector with 
a single-sided SSD (SSSSD) for AF as layer 1, a double- 
sided SSD (DSSSD) for A E; as layer 2 and a 3 x 3 CsI (TI) 
array for residual energy measurements as layer 3. Both of the 
SSSSD and the DSSSD are the BB7 type (2 mm strip width, 
32 strips each side) from MICRON Company. Each CsI (TI) 
crystal is a square pyramid with a dimension of 23 x 23 mm? 
of the front side, 27 x 27 mm? of the rear side and 50 mm of 
the height. The Photo Diode (HAMAMATSU S3204) is used 
to read out the signal from CsI. Fig. 1 presents the detector 
setup of CSHINE in the experiment. For the details and the 
performance of CSHINE, one can refer to [58, 59]. Table 1 
presents the distance d from the center of each detector to the 
target, the polar angle 0, the azimuthal angle ¢ and the sen- 
sitive area S of the SSDTs and the PPACs in the experiment. 
o» The thicknesses of the AE; and A E» for each SSDT are also 
100 listed. 


81 


Si(Au) Telescopes 


PPACI 


Fig. 1. The schematic view of CSHINE Phase-II. 


TABLE 1. Geometric parameters of PPACs and SSDTs. 


Detector d(mm) 6(°) (°) Simm?) AE; (um) AE» (um) 
SSDT1 315.5 18 302 64x64 304 1010 
SSDT2 275.5 25 218 64x64 305 1008 
SSDT3 275.5 31 126 64x64 110 526 
SSDT4 215.5 51 81 64x64 70 306 
PPACI 427.5 50 0 240x280 — — 
PPAC2 427.5 55 180 240x280 — — 


PPAC3 100 180 240x280 


:: The fission fragment detector, PPAC, is a kind of multi- 


wire chamber working in the region of limited proportion- 
ality. The signals induced by the incident fragments on an 
individual wire of the anode plane, either X or Y, are trans- 
ferred through a delay line to both ends. The time delay of 
the two signals X; and Xə (Yı and Y2) with respect to the 
signal collected on the cathode plane, which delivers the tim- 
ing information, gives the X (Y) position of the hit in the 
sensitive area. Fig. 2 presents the schematic view of the me- 
chanics of the PPAC. The total thickness of the sensitive gas 
volume is about 2 cm. The PPACs were operated with 4.5 
mbar isobutane at 465 V voltage. Under this condition, the 
fission fragments can be recorded with efficiency above 95% 
114 but the LCPs and IMFs are suppressed. 

Fig. 3 (b) shows a two-dimensional histogram of Y, — Y2 
us X, — Xə for PPACI as an example. The projections to X 
and Y direction are plotted in Fig. 3 (a) and (c), respectively. 
118 A good performance in timing, corresponding to good posi- 
tion resolution, manifests itself in the sharp boundary for the 
two-dimensional distribution and the well-separated individ- 
ual peaks on the projections. The distance of the neighboring 
wires is 4 mm, and there are 61 peaks and 71 peaks in Fig. 3 
123 (a) and (c), respectively. The time resolution of or = 300 ps 
124 and the position resolution of oy = 1.35 mm can be derived 
125 from the data. The overall performance of the PPACs can be 
12 found in [62]. 
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Fig. 2. (Color online) The schematic view of PPAC as the fission 
fragment detector. 


(a) 


—100 100 0 2000 4000 6000 68000 


0 
X,- X, (mm) Counts 
Fig. 3. (Color online) The scattering plot of Y; — Y2 vs. X1 — X2 of 
PPACI. Panels (a) and (c) are the projections of the two dimensional 
plot of X-Y in panel (b). 


SSDTs are used to measure the LCPs and IMFs in coinci- 
128 dence with the fission fragments. To reduce the total num- 
129 ber of electronics, each neighboring strips are merged in one 
channel, correspondingly the granularity is reduced. Multi 
tracks may fire each SSDT. In order to reconstruct the tracks 
in the SSDTs, a novel algorithm has been developed with spe- 
cial care on the charging sharing effect. More than 80% of the 
hits in all layers of the SSDTs can be recognized and assigned 
to certain tracks. For details, one can refer to [60]. 

The trigger system of CSHINE was designed for both beam 
experiment and calibration. The timing signals of PPACs 
were discriminated by CF8000 and logically calculated by a 
CO4020 to generate the PPAC inclusive signals and PPAC 
two body coincidence signals. The logic hit signals of SSDTs 
were extracted by the front-side of the DSSSD (A E5) with 
the discrimination of MSCF-16, which generates an analog 
Multi-Trig signal proportional to the number of fired strips in 
the same module (16 channels). Both inclusive and exclusive 
logic signals can be generated by discriminating the Multi- 
Trig signal at different threshold settings. In beam experi- 
ment, the trigger signal contains SSD two body events, PPAC 
two body events and the coincidence of PPAC two body with 
SSD one body events. Besides, the inclusive trigger for every 
individual detector were also constructed, optionally turned 
on for the detector calibration before or after the beam data- 
taking. For details, one can refer to [60]. 
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IIl. RECONSTRUCTION OF THE VELOCITY 


153 


We concentrate on the reconstruction of the fission events. 
The flying path of the fission fragments (FFs) can be well 
determined by the PPACs delivering good position informa- 
tion. The velocity of each FF, on the other hand, is derived 
from the timing information. In our experiment, the starting 
time information are provided by the radio frequency (RF) of 
the accelerator. The RF signal, usually a signal in sinusoidal 
form, is discriminated by CF8000 module and input to the 
time-digital-converter (TDC). Generally, for a particle firing 
a given detector, the time of flight (TOF) is written as 
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164 TOF = tact — tre — Cae (1) 

where tact and tgr are the time signals of the detector and 
the RF, recorded by the corresponding TDC channels, respec- 
tively. The unit of both signals are calibrated in ns using pre- 
cise time calibrator. Cae. is a constant representing a fixed 


delay in electronics. The velocity is then computed by 
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170 
where L is the length of the flight path from the target to 
the hit position where the particle fires on the detector. 
In order to verify the validity of the above method to mea- 
sure the TOF, we use the calibrated o particles, where the 
7s velocity is alternatively derived from the energy measured in 
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176 SSDT3. Fig. 4 (a) presents the correlation between the a en- 
17 ergy and the TOF derived by eq. (1). The theoretical curve 
ie fitting the E-TOF profile applies a constant Cact = 431.8 ns. 
179 Panel (b) presents the difference between the TOF measured 
1: by TDC and the value L/v( E), with v being derived from 
181 the energy Ea and L is the distance from the target to hit 
182 position in SSDT3. The width of 1.3 ns is obtained by the 
183 Gaussian fitting. Subtracting the contribution of the energy 
184 uncertainty, the resolution of about 1.0 ns of the TOF is ob- 
185 tained for SSDT. For the FF measured in PPACs, the TOF 
186 resolution will be comparable because the timing resolution 
187 Of PPACS is 300 ps, slightly better than the SSDT. 
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Fig. 4. (a) Correlation of the TOF signal and the energy of o parti- 
cles recorded by SSDT3, and (b) The distribution of the time differ- 
ence between the TDC measurement and the calculated value from 
the calibrated energy. 


With the validity eq. (1), we can discuss the TOF of the 
FFs. Since PPACs can not identify the charge or mass of the 
FFs, nor the total kinetic energy, one relies on the determina- 
tion of the velocities which requires the TOF information and 
the lengths of the flight path L. In eq. (1), the relative differ- 
ence of the delay constants, Cppaci — Cppac2, between two 
PPACs can be adjusted to zero using pulser prior to the exper- 
iment and the systematic uncertainty can be well controlled 
within +2 ns. However, the absolute value of Cppac of each 
individual PPAC can not be determined like in the SSDT, be- 
cause the particle type and the total energy is unknown. To 
overcome this difficulty, we use the Viola systematics that the 
relative velocity of the FFs is averagely 2.4 cm/ns [63] . And 
201 hence, by tuning the constants Cppac, one can optimized the 
202 value at (vpp) = 2.4 cm/ns. Fig. 5 presents the distribu- 
tion of upp at different delay constants. One can readily see 
that varying the delay constant by 1 ns, the peak position of 
vpr moves significantly. In our experiment, Cppac = 115.5 
ns is optimized, and the corresponding distribution of upp is 
207 plotted in panel (b) . Known from transport model calcula- 
tions, the variation of (vrp) is better than 0.1 cm/ns, a sys- 
tematic uncertainty of 2 ns of the TOF of FFs is estimated, 
on which the conclusion of the following analysis will not 
21 Changed. We recall that the coincident events recorded by 
22 PPACI with PPAC2 (marked by PPACI x 2) or with PPAC3 
23 (marked by PPACI x 3) are the fission fragments , since the 
214 HV condition is set that the response of PPAC to energetic 
2:5 LCPs and IMFs are totally suppressed, given that the energy 
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loss of these particles is lower than FFs by more than one 
order of magnitude. In addition, the correlation between a 
heavy projectile-like fragment (PLF) and a target-like frag- 
ment (TLF) is out of the current geometrical coverage. 
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Fig. 5. The relative velocity distribution vfr of the two fission frag- 
ments with different parameter setting of Cppac. 


IV. RESULTS AND DISCUSSIONS 


Before discussing the reconstruction of the fission events, 
we first define the kinetics of the fission. In the picture of in- 
complete fusion, a heavy TLF is formed in the fusion of part 
of the projectile and the target. The fraction of the momentum 
of the projectile transferred to the TLF is called linear mo- 
mentum transfer (LMT). With a certain probability depend- 
ing on the total angular momentum of the reaction system, the 
TLF may undergo fission or fast fission in competition with 
the emission residue channel. For the fission events, Fig. 6 
presents the kinetic geometry of the TLF fission event. As 
shown, the origin point O is the target nucleus in the labo- 
ratory system, the vector OO represents the direction of the 
beam. The velocity vectors vr, and v, of the two fission frag- 
ments in the laboratory system are represented by OA and 
OB. vy is the velocity of the TLF, and the velocities of the 
two fragments in the center-of-mass system of the fission- 
ing TLF are represented by Uy and Te, sitting back-to-back 
collinearly. Here, we define plane OAB as the fission plane, 
plane OO A as the projection plane, and the reaction plane is 
defined as plane DOO’. 

Once the TOF is determined, the velocity of the FFs can 
be computed event-by-event using the hit positions of the FFs 
in PPACs. Then the whole fission event can be reconstructed. 
Fig. 7 presents the distribution of the velocity of the fission 
fragments in PPAC1 x 2 and PPAC1 x 3 events, respec- 
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Fig. 6. The geometric diagram of the velocities of the two fission 
fragments from TLF. 


tively. Here the vr, is for the FF recorded in PPAC1 and vp, is 
for the FF in PPAC2 (PPAC3) in PPACI x 2 (PPACI x 3) 
events. It is shown that for the PPACI x 2 events, since the 
two PPACS are nearly symmetric w.r.t the beam, the distribu- 
tions of vr, and vr, are very similar. Meanwhile, high velocity 
tail is evidently presented in the PPACI x 2 events and the 
velocity spectra are wider than that in PPACI x 3 events. 
This component is mainly due to the events with smaller fold- 
ing angle corresponding to larger LMT, as will be discussed 
below. 
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Fig. 7. The velocity distribution of the two fission fragments in 
PPACI x 2 and PPACI x 3 events. 


After the velocities of the two FFs are determined, the fold- 
ing angle method can be applied to calculate the LMT of the 
reaction. Recalling the definitions in Fig. 6, the folding angle 
Orr is defined as the angle ZAOC, spanned by the projec- 
tion of the velocity vectors of the FFs on the projection plane. 
Clearly, given the velocities vr, and vr,, the folding angle de- 
pends on the velocity of the TLF vy, i.e., the larger vq, the 
smaller Opp. Fig. 8 presents the distribution of the folding 
angle Opp. The coincident events of PPACI x 2 distribute 
in the range between 70? — 120? with the peak situating at 
about 95?, corresponding to a larger LMT (red), while the 
PPACI x 3 events sit in the range of 120? — 170? corre- 
sponds to a smaller LMT (blue). The valley between the two 
components is simply due to the deficiency caused by the gap 
between PPAC2 and PPAC3, and the efficiency arising from 
the incomplete azimuthal coverage is not corrected in the plot. 
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Fig. 8. (Color online) The folding angle distribution. 


As the fission geometry is determined by the two velocity 
vectors, it is of interest to look at the planarity of the fission 
events. Fig. 9 (a) present the azimuth correlation of the two 
FFs by scattering plot between the azimuthal angle difference 
A and the folding angle Orr. It is clearly shown that for the 
FFs from both the central and peripheral reactions, the most 
probable value situates at Ao = 180° following the picture 
that the system undergoes a binary decay. Here we note that 
the Aó is a directly measurable quantity relying on no as- 
sumption. The inset displays the projection distribution of 
A4, a standard deviation of about o(A¢) z 10? is derived. 
Such broadening suggests that the emission of LCPs or IMFs 
may change the flight direction of the FF and smear the back- 
to-back feature of the fission event. In order to see its evolu- 
tion with the violence of the reaction, we plot in panel (b) the 
standard deviation of azimuthal angle c(Aó) as a function 
of folding angle Opp. It is shown that o (A4) decreases with 
Orr in the whole range , with an exception near Opp ~ 130°, 
where a discontinuity appears due to the gap between PPAC2 
and PPAC3. To exclude the possible reason that this trend 
originates from the asymmetry of the geometrical locations of 
the PPACS, we restrict further the analysis on the events with 
the two FFs flying symmetrically to the beam direction, i.e., 
with the condition of 0; = 62, where 0; is the polar angle of 
the it! fragment. The result is depicted by the black squares, 
where the bin width of each Opp is +2.5°. This condition 
is applicable only for the PPAC1 x 2 fission events because 
these two PPACs are placed in approximate left-right sym- 
metry with respect to the beam line. Clearly seen, the data 
points with the symmetry condition are sitting on top of those 
without the condition, suggesting that the decreasing trend of 
c (Ado) as function of Opp is truly due to the reaction vio- 
lence. Since the post-scission particle emission changes the 
velocity of the FF due to recoil effect, the trend suggests that 
in the fission following the intermediate energy heavy ion re- 
actions, there is sufficient excitation energy left at the scission 
point depending on LMT. In the reactions with larger LMT, 
more excitation energy is left and released through the parti- 
cle emission in the post-scission stage. It is consistent with 
the picture of fast fission, instead of statistic fission in which 
the excitation energy is nearly depleted at the scission point. 


The dynamic feature of the fast fission can be further ex- 
plored from the velocity distribution of the FFs. We present 
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Fig. 9. (Color online) (a) The azimuthal angle difference Ad as a 
function of the folding angle Orr of the fission fragments. The inset 
displays the projection distribution of Ao. (b) c(Ad) as a function 
of Orr ^ 


in Fig. 10 the average (a) and the standard deviation (b) of 
the velocities of the FFs recorded in the PPACS as a func- 
tion of folding angle. From panel (a), it can be seen that 
the average velocity value (vr) decreases with the folding an- 
gle. Panel (b) presents the standard deviation of the velocity 
C (vy), and it is also clear that the broadening of the velocity 
of the FFs decreases with the folding angle. This result is con- 
sistent with the picture seen in the trend of c(A4) in figure 9. 
The scission point is early reached when the excitation energy 
of the fissioning TLF is still high, thus, the statistical fluctu- 
ation (corresponding to the left excitation energy) enhances 
the variance of the velocity of the FF. This is consistent with 
the earlier experimental observation in Ar +7°9 Bi reactions 
at 25 MeV/u [64, 65]. 

Finally, the dynamic feature of the fast fission may also 
cause the anisotropy of the angular distribution of the fission 
axis. Here the fission axis is defined as the vector of the rel- 
ative velocity Upp from the fọ to fj, where f, and f; are the 
two fragments. The angle of the fission axis Arr is the an- 
gle of the fission axis w.r.t. the beam axis, as defined in [55]. 
Usually, the experimental measurement of the distribution of 
Apr requires fine correction of the geometric efficiency, and 
hence it is better feasible using 47 detecting system. 

In our experiment, PPACs cover only part of the whole 
space, in order to introduce less ambiguity to the geometry 
efficiency correction, we fix the direction of the first FF in 
PPACI. In this case, we need only to correct the efficiency of 
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Fig. 10. (Color online) The average value of velocity (a) and its 
standard deviation (b) of the velocity distribution as a function of 
folding angle. 


342 the second FF on PPAC 2 and PPAC3, and the trend of the an- 
343 gular of the fission axis can be inferred. Figure 11 presents the 
344 correlation plot of parallel and transverse velocity of the FFs. 
345 Here the transverse velocity of the FFs recorded in PPACI is 
» s defined as positive and that in other two PPACs are defined as 
347 negative. It is clearly shown there is a dead area of less than 
348 20? between PPAC2 and PPAC3. The dashed lines define a 
34 narrow range of 40° < 0; < 45°, which is fixed in the in- 
sso Vestigation of the distribution properties of the fission axis as 
«351 following. 
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Fig. 11. (Color online) The transverse and parallel velocity distribu- 
tion of the fission fragments. 


The distribution of do /d cos(Arr) is plotted in Fig. 12. 


352 


sss The geometric efficiency arising from the incomplete azimuth | se 


354 coverage of PPAC2 and PPAC 3 is corrected in each Apr 


3s bin. The events of PPACI x 2 and PPACI x3 are rep- 
sss resented by symbols, while the curve represents the sum. 
357 It is clearly shown that the deficiency of the gap between 
358 PPAC 2 and PPAC 3 causes a kink in a wide range in 63° < 
39 App < 80°. Regardless of the kink area and the uncovered 
3 region within App < 50°, it is shown that the distribution of 
361 dN/dcos(App) increases steadily with cos(App) and tends 
32 to peak at forward angle, which is at variance with the expec- 
33 tation of isotropic distribution for statistic fission. The trend 
34 is in qualitative agreement with the previously reported ex- 
365 perimental results in HIR at Fermi energies [48, 49]. 
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Fig. 12. (Color online) The angular distribution dN/dcos(Arr) of 
the fission axis w.r.t. the beam. 


V. SUMMARY 


366 


In summary, the fission fragments in 25 MeV/u 99Kr 
-?95pb reactions have been measured with CSHINE detec- 
tion system. At current phase three PPACs and four SSDTs 
are mounted to measure the fission fragments and the coinci- 
dent LCPs and IMFs, respectively. Using the timing and posi- 
tion signals of the PPACs and the start timing from the RF of 
the accelerator, one can measure the velocities of the fission 
fragments and reconstruct the fission events, where the lin- 
ear momentum transfer can be derived by the folding angle. 
It is shown that, the width of the azimuthal angle difference, 
the mean value and the width of the velocity distribution of 
the fission fragments decrease with the folding angle. The 
anisotropic angular distribution of fission axis has been ob- 
served. The results are consistent with the picture that fast 
fission occurs. Prospectively, with the ability of reconstruct- 
ing the fission events, CSHINE provides the opportunities in 
the studies of isospin dynamics and nuclear symmetry energy 
by further counting the coincident isotope-resolved LCPs and 
IMFs. 
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